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Abstract. In this paper we study parameter families of truncated horseshoes 
as models of multiscattering systems which show a transition to chaos with- 
out losing hyperbolicity, so that the topological features of the transition are 
completely describable by a parameterized family of symbolic dynamics. At 
a fixed parameter value the corresponding horseshoe represents the set of or- 
bits trapped in the scattering region. The bifurcations are a pure boundary 
effect and no other bifurcations such as saddle center bifurcations occur in this 
transition scenario. Truncated horseshoes actually arise in concrete poten- 
tial scattering under suitable conditions. It is shown that a simple scattering 
model introduced earlier can realize this scenario in a certain parameter range 
(the "truncated sawshoe") . For this purpose, we solve the inverse scatter- 
ing problem of finding the central potential associated to the sawshoe model. 
Furthermore, we review classification schemes for the transition to chaos of 
truncated horseshoes originating from symbolic dynamics and formal language 
theory and apply them to the truncated double horseshoe and the truncated 
sawshoe. 
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1. Introduction 

We are going to study the following type of scattering systems. A set of scat- 
terers is arranged in some region of the real plane according to some deterministic 
law. A test particle is injected into this scattering region and will undergo multi- 
ple scatterings before (if ever) it will escape. A discretization of this dynamics is 
assumed to be possible in the form of a map T, which describes a single scatter- 
ing event. Certain features of this scattering system such as the singularities of 
reaction functions or cross sections, can be described by the set of trapped tra- 
jectories which never escape from the scattering region. Introducing and varying 
an external parameter, controlling the potential strengths say, varies the trapped 
set of the scattering system. In 2 D the trapped set typically takes the form of a 
Smale's horseshoe. 

In this paper we suggest the following abstraction. We interpret parameter 
families of truncated horseshoes as models of multiscattering systems which show 
a transition to chaos completely describable by a parameterized family of sym- 
bolic dynamics 1 (SD). At a fixed parameter value the corresponding horseshoe 
represents the trapped set. Deforming or chopping the horseshoe will change the 
trapped set and consequently the physical properties of the whole system. 

A horseshoe can be generated or destroyed in two ways: either by chopping 
off parts of it or by smoothly deforming and disentangling the horseshoe tangle 
yielding truncated horseshoes and disentangled horseshoes, respectively, (cf. fig- 
ure 1). The property distinguishing these two scenarios is that truncation of a 
horseshoe preserves the hyperbolicity 2 of the underlying dynamical map, whereas 
smooth deformations produce tangent bifurcations with the consequent loss of hy- 
perbolicity. This is why truncated horseshoes are the simpler models, which lend 
themselves more readily to a rigorous approach. 

Disentangled horseshoes as models for the transition to chaos of potential scat- 
tering have been discussed by Bleher, Ding, Grebogi, Ott, and Yorke in a se- 
ries of papers (cf. [2, 1, 6]). But truncated horseshoes also arise in concrete 
potential scattering, if one chooses a suitable arrangement of non-overlapping, 
attractive potentials with compact supports (i.e. of finite range) and permits 
non-differentiability of the potentials at the boundary of their supports. A sim- 



1 Some basic concepts about SD can be found in appendix A 

2 A set A is called hyperbolic for a map or semi-flow T if it is compact, invariant under T and 
its tangent bundle can be continuously decomposed into a direct sum of stable and unstable 
manifolds, whose contraction factors have a uniform bound < 1. 
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FIGURE 1. (a) The time 2 and -2 image of a complete Smale horse- 
shoe together with a coding for the vertical (forward) and hori- 
zontal (backward) strips, (b) A smooth deformation yielding an 
incomplete Smale horseshoe, where not all forward strips intersect 
all backward strips, (c) A Smale horseshoe truncated by a family 
of level lines of a pruning function /. The shaded convex region is 
supposed to satisfy f(u) > n. The same 4 second generation rect- 
angles are lost as in Fig. 1 (b). The forbidden symbol combinations 
are {0001,1001,1011,0011}. 
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pie example of this kind (cf. [15]— [17] ) consists of such potentials which are 
additionally identical and central with a (single) sawtooth as deflection function 
and which are arranged at unit distances along a straight line. The potentials 
are parametrized by the slope k > of the deflection function and by the range 
parameter < 6 max < 0.5. 

We shall show in the present paper that a sawtooth deflection function in- 
deed arises from a suitable attractive central potential. The jump discontinuity 
of the sawtooth proves to have no severer effects than different right and left 
derivatives of the potential at the boundary of its support. Everywhere else the 
potential is smooth. Moreover, the trajectories are everywhere once continuously 
differentiable. 

Actually, this model also has a non-hyperbolic chaotic regime for < k < 4, 
where the horseshoe develops by disentangling (cf. [11]). What we are inter- 
ested in here however, is the parameter range for which a transition to chaos 
occurs within the hyperbolic regime (6 max small enough and k > 4). One must 
be aware that hyperbolicity alone does not assure structural stability and thus 
does not exclude bifurcations. In the systems we are concerned with here, there 
are boundary effects originating from the non-differentiability of the attractive 
potentials on the boundary of their support. These boundary effects generate 
the bifurcations leading to a transition of the dynamics from regular to chaotic 
behaviour (see next section). Only for the purpose of analysing the transition 
more precisely, we shall introduce a linearization. The qualitative features of 
the linearized version and the original model are the same. As the horseshoe 
associated to this model resembles a truncated sawtooth in the asymptotic limit 
and in the linearized model, we speak of a "sawshoe" in order to distinguish this 
particular model. 

The advantage of truncated horseshoes is that their transition to chaos reduces 
to the problem of analysing the evolution of the trapped set in terms of the 
evolution of an associated SD. 

In the family of associated SD the evolution of the trapped set occurs as an 
evolution of the set of "permitted" symbol sequences each describing an orbit of 
the trapped set. This set of permitted symbol sequences is called a subshift. If it 
comprises all symbol sequences it is called a full shift. As one often represents a 
subshift as a tree of permitted sequence segments, the evolution of subshifts cor- 
responds to an evolution of trees called pruning process. Multiscattering families 
which are described by a pruning process will be called pruned systems. 

We have identified two scenarios by which pruned systems can become chaotic: 
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a transition via a cascade of bifurcations and an abrupt transition. 

Both, chaoticity and topological entropy of a pruned system essentially mea- 
sure its size (cf. definition 3.4 and remarks there). Not even topologically do they 
determine the pruned system completely. In symbolic dynamics, there is a fur- 
ther refinement of classification by the Markov property of permitted sequences, 
namely into finite and infinite subshifts. 

A means to classify infinite subshifts further is provided by formal language 
theory. A formal language is just a set of finite strings (as opposed to sequences 
in the case of subshifts) of symbols taken from a finite set called the alphabet. 
By relating formal languages with symbolic dynamics one can use formal lan- 
guage theory to study the structural complexity of pruned systems. Structural 
complexity means roughly the type of rules needed to describe the production 
of those words belonging to the formal language in question. The theory of au- 
tomata and formal languages shows that there are types of rules which are not 
arbitrary but have certain universality properties. The best known examples are 
the Turing machines and the finite state automata, both of which belong to the 
Chomsky hierarchy of automata (or equivalently of languages). Important now 
for the applicability to dynamical systems is that the Chomsky hierarchy offers 
a classification scheme, which proves to be a also a topological invariant, i.e. 
topologically equivalent SD fall in the same complexity class. 

There is now a bridge between structural complexity and the onset of chaos in 
the following sense: The onset of chaos in pruned systems shows certain features 
of a phase transition, seen for instance in the topological entropy as a function 
of the family parameters (cf. [16]). Recently there have been a few studies 
about the complexity theoretic description of phase transitions based upon the 
informational diversity and computational complexity of observed data (cf. [4]). 
It was suggested there that it is at phase transitions that high level 3 computation 
occurs, so that computational ideas provide a new set of tools for investigating 
the physics of phase transitions. 

1.1. Overview. The contents of this paper is as follows: Section 2 explains, 
why truncated horseshoes can model chaotic scattering. Two simple examples 
are treated, the double horseshoe and the sawshoe. This part is based on [15] 
and [17]. A continuous and but for the boundary of its support smooth potential 
associated to the sawshoe model is given. This is new. The details are worked 
out in appendix B. 



3 e.g. in the sense of the Chomsky hierarchy. 
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Sections 3 and 4 review parts from [17], [19] and [10] for section 4.1. Section 3 
introduces formal multiscattering systems (FMS) (definition 3.1) and quotes suffi- 
cient conditions for them to have a symbolic dynamics (theorem 3.1). Truncation 
of FMS is then determined by a set of level lines of a cut-off function. Whenever 
a particle overjumps such a level line it is supposed to escape from the scatter- 
ing region (definition 3.3). Furthermore, two scenarios for transitions to chaos of 
truncated horseshoes are discussed in this section: the cascade transition with in- 
finitely many bifurcations preceding and accumulating at the critical point where 
chaos sets in (this transition b strong resemblance to the period-doubling 

scenario of bounded dynamical systems) and the abrupt transition where "1 + e" 
bifurcations suffice to generate chaos, i.e. the dynamics restricted to the trapped 
set is not yet sensitively dependent on initial conditions at the critical value, but 
at any e beyond it is chaotic. These concepts are applied to the double horse- 
shoe and the sawshoe, which are shown to follow the cascade and the abrupt 
transition, respectively (corollary 3 and theorem 3.2). 

Section 4 discusses the same transitions from the point of view of complex- 
ity theory. To this purpose the notion of dynamically generated languages is 
introduced. These are languages which can be naturally associated to subshifts 
(proposition 4). The theory of formal languages offers the concepts to grasp 
the observation that typically the rules specifying the permitted sequences in a 
pruning process become more and more complicated close to the critical point of 
a dynamical system, indicating the need of a higher level organization of these 
rules. The sine qua non for the application of theses ideas from formal language 
theory is given in theorem 4.1. This theorem essentially states that the com- 
plexity measure most often used in formal language theory, namely the Chomsky 
hierarchy, has an intrinsic dynamical meaning, i.e. is not an artifact of the special 
description, which means here the chosen symbolic dynamics. The complexity 
transitions of the two example horseshoes are described by corollary 4 and theo- 
rem 4.3. 

Appendix A and C contain some basic concepts from symbolic dynamics and 
automata theory, respectively; in appendix B a potential with sawtooth deflec- 
tion function is calculated and some of its properties, such as uniqueness, are 
discussed. A Method from classical inverse scattering theory (cf. [12]) is used. 

1.2. Acknowledgement. I would like to thank Uzy Smilansky and Andreas 
Knauf for discussions. 
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2. Truncated horseshoes as scattering models 

2.1. The truncated double horseshoe. One of the first models studied in 
the days of the rediscovery of chaotic scattering by physicists in the eighties (cf. 
[7]), was the three disk model: A particle is scattered in 2 D at three hard disks 
arranged at the vertices of an equilateral triangle. When one injects a whole 
bunch of particles coming from an interval which together with a given direction 
of motion form a strip enveloping one of the disks, one notices that those particles 
staying in the scattering region of the three disks for another scattering incident 
originate from two smaller subintervals. These subintervals in turn contain two 
still smaller subintervals each, which represent particles realizing at least one 
more, that is at least three, scatterings. Iteration of this procedure generates the 
Cantor set of initial conditions, for which particles will be trapped forever in the 
three disk system. 

In potential scattering the more natural space to study is phase space. We 
restrict our attention to 2 D, The canonical dynamics associated to a Cantor set 
in 2D is Smale's horseshoe. Now, in order to find the most simple multiscattering 
system, we are going to turn the problem upside down: 

Is there a (natural ?) way to interpret the abstract Smale horseshoe as a multi- 
scattering dynamics ? 

I suggest the following interpretation: Take the unit square Q\ as representing 
the part of phase space which has the right properties to be scattered at least 
once at a single potential of some hypothetical arrangement of potentials. We 
shall call this part of phase space here fundamental domain. Let the dynamics 
of this single scattering be given by Smale's horseshoe map Tj 5 , which stretches 
and folds the unit square back onto itself. Those parts of the square's first image 
that are mapped outside the square (for example the bent part of the horseshoe), 
is supposed to represent the phase space points of particles which have escaped 
the scattering region for good. Whereas the two strips within the unit square 
represent particles which will achieve at least one more scattering. Iteration yields 
the forward trapped set of all points which have undergone an infinite number 
of scatterings. Backward images yield in the same way the backward trapped 
set of all points which will undergo infinitely many scatterings. The Cantor set 
of Smale's horseshoe is then just the intersection of the forward and backward 
trapped set, called the trapped set KiT^Qi). By labeling the two strips of the 
horseshoe, its Cantor set can be identified with the set of all sequences with 2 
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symbols, say and 1, the so-called full shift of two symbols. 4 

The arguably simplest way to generate more complicated trapped sets than 
the full shift of two symbols, is by truncating the unit square horizontally at a 
point < y max < 1- The new fundamental domain is therefore: 

(1) Q, max = {(^,y)G[0,l] 2 ;0<y<y max } 

Points mapped outside Q Vra ^ x are supposed to have escaped from the scattering 
region. The resulting map is denoted by T^ max . The lowering of y max below 1 has 
the effect that not all sequences of two symbols correspond to trapped orbits : this 
is called pruning. The question which will occupy us throughout this paper, reads 
as applied to the truncated horseshoe: How can the trapped set A(T^ max , Q yraa _ x ) 
be characterized as a function of the family parameter y max . 

The truncated Smale's horseshoe is not yet the most simple imaginable mul- 
tiscattering system, because it contains a reflection. 5 We therefore take instead 
(see figure 2) Smale's double folded horseshoe map T 2/max (which as it happens was 
the original form found by Smale). The double horseshoe family with trapped 
set 

oo 
i— — oo 

= {u = (x,y) G Qi : ir y (Tl(u)) < y max G Z} 

where tt v is the projection onto the y-component, is equivalent to the symbolic 
dynamics (cr, E^) over the alphabet A = {0,1} with trapped subshift E„ and 
truncation parameter fi gA No , determined by the fundamental domain: 

(2) 1^ = {a gA z : P(a)</i} 
and the truncation (pruning) function 

oo 

(3) P{a) = J2^2- 1 

8 = 

For an orbit a gA z /<t we define 

(4) P(a) := sup P(a) 

aGa 



4 Some pertaining basic definitions of finite, infinite, cyclic subshifts and (semi-) conjugacies 
can be found in the appendix. 

5 The eigenvalues are negative in one strip and positive in the other. 
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FIGURE 2. (a) The twice folded, piecewise linear horseshoe map 
applied to the unit square Q = Q\. The symbolic dynamics chosen 
uses the first generation horizontal strips. (b)The horizontally trun- 
cated double horseshoe. As the reflection has been removed, the 4 
second generation rectangles lost for the chosen cut-off parameter 
i/max correspond to the symbol combinations 0011, 0111, 1011, 1111, 
respectively. 
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Observe that we interpret sequences with nonnegative indices such as ji as binary 
numbers, i.e. in the base 2. 



2.2. The truncated "sawshoe" as a potential scattering model with 
complete SD:. We address in this section a multiscattering model introduced 
by the author and U. Smilansky in [15]. Its salient feature is that it yields 
a parameterized family of discrete dynamical systems which remain expansive 
throughout a transition to chaos. Therefore this transition to chaos can be de- 
scribed at least topologically in a complete way by a family of symbolic dynamics 
(cf. theorem 3.1 of section 3). 

Suppose a particle of a fixed asymptotic energy E > is scattered at an 
attractive central potential V = < 0. For the sake of simplicity we assume 
V to be of finite range 6 max > 0. Conservation of angular momentum restricts 
the dynamics to a plane in configuration space. The equations of motion can 
be integrated to yield the polar angle (j> of the trajectory as a function of the 
impact parameter b and the radial distance r. The asymptotic final angle 

(for t — > oo) is expressible as a function of the asymptotic initial angle x^ (for 
t — y — oo) and the impact parameter as: 

(5) xf = xf ] + 9 v (b) 

where the deflection function 6y is given by equation (31). Multiscattering can 
then be achieved by positioning in configuration space copies of V according 
to some group action and avoiding overlap. The particle ejected from may 
either escape the scattering region for good or it may approach the next potential 
V( 2 ) under a new initial angle and a new impact parameter both of which are 
functions given by the geometry of the arrangement of the potentials. The map 
(xf'\b^) i — y (x^^b^) is then the discretized scattering map. 

Placing the potentials along a line at equidistant positions yields the following 
scattering map: 

T : S 1 xR — y S 1 xR 

(6) (x\ ( x + e v {b) \ 

\b)^ {b+G(x + e v (b))J 
where S 1 is the unit circle and 

(7) G(x) = — sgn(cos x) sin x 
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FIGURE 3. The potential Vz K , which has the sawtooth as de- 
flection function for some 6 max - and A;- values. 

In order to arrive at a concrete model we choose as deflection function By the 
sawtooth function k > 0: 



In the papers quoted the question whether there is indeed a potential with such 
a deflection function was left open. In the appendix this question is answered 
positively. It is shown that is the deflection function of a potential Vz K which 
is smooth in ]0, b max [ and continuous on IRq . The potential Vz K is defined by the 
following parameter representation: 



together with Vz K (r) := if r > 6 max - Figure 3 shows this potential for some 
values of k and 6 max - Some properties are discussed in appendix B. Figures 4 
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FIGURE 4. Some sample trajectories in the potential Vz K before 
the background of a density plot of the potential. In the right figure 
the trajectory winds 3 times around the center. 

(a,b) show some sample trajectories in the potential Vz K - 

Thus, the family Tkfi ma , x together with a parameter set V C Rx]0,0.5[ defines 
a 2-parameter family of potential scattering systems. The associated family of 
trapped sets is 

= {u = (x,b) g <5 6max ; K P (r^ 6max (u))| < 6 max VieZ} 

= {u G Afc,6^ ax ; \n (tI(u)) I < LxVi G Z} for any b' max > 6 max 

where the pruning map ttj, is the projection onto the 6-component. Varying 
only the parameter k will lead to a deformation of a given trapped trajectory, 
which then may miss one of the unchanged scattering potentials and become a 
free trajectory (cf. figure 5 (a)). Varying only 6 max changes the selection rule 
\b\ < fe max , so that a formerly trapped trajectory might miss one of the shrunken 
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FIGURE 5. Example of (a) the influence of a variation of the pa- 
rameter 6 max on a trajectory with fixed impact parameter b and (b) 
of a variation of the parameter k. 
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potentials and again become free (cf. figure 5 (b)). In both cases the trapped 
set changes. The difference between the two cases is that varying only 6 max does 
not change the scattering map restricted to the shrunken interaction set, i.e. for 
i>L a ^ < &JL V we have T kh 2 \Q h i = T kh i \Q h i . Figure 6 shows the truncated 

max max ^i u max 1 ^ u max ^i"max 1 ^ u max ° 

sawshoe. 

Now we shall introduce another simplification. From now on we restrict our- 
selves to the piecewise linearization of the scattering map given by 

\ —x if \x\ < it/2 

^ ^ = 4. -f I I Icy 

I —X + 7T II \X — 7T| < 7T/Z 

This linearization corresponds to the limit 6 max — > in the physical scattering 
family and describes therefore its asymptotic behaviour. In the parameter set V 
defined as 

(11) V = {(k, 6 max ) ; < 6 max < 0.5, k > 4, tt < k 6 max < 2ir} 

we observe a transition from a regular to a chaotic dynamics via a pruning process, 
as will be made clearer presently. 

A (generalized) symbolic dynamics for the sawtooth family which yields a sim- 
ple symbolic form of the pruning map ttj, is the following: 

The sawtooth family is 2:1 semi-conjugate to the following pruned family over 
the alphabet A = {-1, 0, 1}: 

(12) with \i G £ N 

whose dynamical map a = —a is the negative shift map (i.e. still a trivial cellular 
automaton). The trapped a -spaces = A(<r, T^) are determined by polynomial 6 
pruning in the alphabet ti = {-2, -1, 0, 1, 2}; 

(13) T, = {aeA z : \P(a)\ < /i} 

P (a) = a ; P t (a) = a_ t + a t} i > 1 

with pruning map 

oo 

(14) P(a) = Pi{a)5~\ P(a) := sup P(a) 



i = 



aGa 



3 It is called polynomial because the pruning functions Pi are all polynomials. 
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FIGURE 6. The "sawshoe" A/ ;j 6 max generated in phase space by the 
scattering system (Tfc,& max , Q& max ) with interaction set Qb rn!ix shown 
in figure (a) and (b) is incomplete for the right figure (b) with the 
smaller 6 max (and consequently smaller Qb ma , x )- 
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Observe that we again interpret sequences with nonegative indices as real num- 
bers, here in the base 5. 

Remark 2.1. The truncation scenario of the sawshoe does not depend on the lin- 
earization (10). The nonlinear dynamical system determined by (7) and (8) is 
hyperbolic over an open parameter subset and undergoes there a hyperbolic tran- 
sition to chaos by truncation of its horseshoe if one restricts to k > 4 and 6 max 
small enough, because hyperbolicity is stable with respect to smooth perturba- 
tions and indeed the sine function converges uniformly toward its linearizations 
around its zeros for 6 max — > 0. One case of this claim is even trivial, namely if 
one fixes a k > 4 big enough and 6 max small enough so that both the system is 
hyperbolic and the symbolic dynamics is the full shift A z and if subsequently one 
starts to decrease 6 max down to zero. Obviously the hyperbolicity is conserved 
under such a parameter variation, because the scattering map changes only for 
the escaping points. Observe that trajectories change their form only within the 
potential's support (see figure 5) and will be pruned if and only if the potential's 
range shrinks below their largest impact parameter. The difference between the 
linearized and the nonlinearized system is that for the latter the pruning map P 
is more complicated due to a loss of symmetry (cf. [18]). 

3. Transition to Chaos of formal multiscattering systems 

In this section we are going to introduce the concept of a formal multiscattering 
system 7 (FMS) in order to isolate the essential features of the examples given 
above. We shall give sufficient conditions for the existence of a SD for a FMS and 
we shall show how a FMS can become chaotic by truncation of the associated 
horseshoe. Finally we shall apply these concepts to the examples introduced 
above. 

3.1. Formal Multiscattering systems (FMS). 

Definition 3.1. Let M be a metrizable space and T : M — > M be a bijective 
map. Furthermore, let Q C M be a compact subset and define the following 
images: 

(15) r := T\Q) n Q, i G Z 

7 In computer science (cf. [13]) a similar concept neglecting however topological structures 
is used and has recently been applied to the study of dynamical systems, namely the so-called 
predicate transformers. 
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The pair (T, Q) is called a formal multiscattering system (FMS) if the following 
conditions are satisfied: 

(i) homeomorphism property: The restriction T\Q — > T(Q) is a homeo- 
morphism 8 . 

(ii) escape property: For all n G N : T n+1 (I- n ) \ Q ^ 0. 

(iii) trapping property: For all n G N : T n+1 (I- n ) nQ^0. 

(iv) no-return property: if x G Q, T(x) ^ Q then Vn G N : T n (x) ^ Q. 

From the definition one concludes that he trapped set A of a FMS is non-empty 
and compact. 

A sufficient property for the existence of a symbolic dynamics (SD) for a FMS 
is expansiveness: 

Definition 3.2. (i) A FMS (T,Q) is called expansive if T|A is expansive in 
the usual sense, i.e. if there is a 8 > Os.t. x ^ y implies the existence of 
a n G Z with d(T n x, T n y) > 8. 
(ii) A finite open cover a of A(T, Q) is called a generator for the FMS (T, Q) 
if V sequences G a z : OT^Ai) contains no more than 1 point in Q. 
If f\T l (Ai) contains no more than 1 point, a is called a weak generator. 

Examples include the truncated double horseshoe and sawshoe from above (be- 
cause they are hyperbolic) and the three disk scattering model because reflection 
at disks defocuses trajectories. 

As in the standard situation of a bounded dynamics (cf. [20]) the existence 
of a generator is tantamount to expansiveness. As one can use a generator to 
construct a SD, expansiveness implies the existence of a SD. More precisely, one 
has the following three statements: 

Lemma 1. The FMS is expansive iff it has a generator iff it has a weak genera- 
tor. 

Theorem 3.1. If (T,Q) is an expansive FMS, then there is a surjective semi- 
conjugacy $ : E — > A(T, Q) from a subshift E over some finite alphabet A. 

Corollary 1. If the FMS (T,Q) has a disjoint generator, then $ can be cho- 
sen to be a topological conjugacy. In this case the trapped set A(T, Q) is totally 
disconnected. 

8 i.e. bijective and both T and T _1 are continuous 
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Finally we are going to make the notion of truncated families of FMS more 
precise: 

Definition 3.3. A truncated family of FMS is given by a parameterized family 
of FMS (T K , Q K )neJ over some interval J C IR n as parameter set, where each Q K 
is the truncation of a common set Q D Q K by a family of level lines given by a 
cut-off function / : M — > IR and an evaluation function e : J — > IR + : 

(16) Q K = {ueQ : \f(u)\<e( K )} 

3.2. Transition to Chaos. Expansiveness of a FMS is not quite enough for 
chaos because A may contain isolated points: 

Definition 3.4. Let (T, Q) be a FMS with trapped set A. We call (T, Q) weakly 
chaotic if its topological entropy 9 H top > 0. The scattering dynamics is called 
strongly chaotic (or chaotic in the sense of Devaney) if (cf. [5]): 

(i) the restricted map T|A is sensitively dependent on initial conditions (s.d.i.c.). 10 , 

(ii) there is a dense orbit in A and 

(iii) the set of periodic points Pery(A) is dense in A: 
LL T (A) := Per T (A) = A. 

In [3] it is shown that part (i) can be replaced by requiring A to have infinitely 
many points. 

We call the opposite of s.d.i.c. stability. If T\A(T,Q) is stable everywhere on 
the trapped set, which means in this context that A is discrete, then we speak of 
a stable or regular dynamics. 

Whereas the symbolic dynamics of a full shift with number of symbols $A > 1 
is always strongly chaotic and that of a finite subshift is strongly chaotic iff the 
transition matrix is irreducible and aperiodic, there is no simple criterion for 
infinite subshifts. 

The following questions about pruned multiscattering families and in particular 
about their symbolic models are of interest for us in this connection: 

• Where is the map ji i— > not locally constant, i.e. where does the 
trapped set acquire new points ? These values will be called bifurcation 
values forming the bifurcation set V = V(r, T). 

9 As A is compact, it is discrete iff it is finite. Hence the topological entropy 
H top (T, Q) := linir-^oo lo s#{" eA ' per T (ti)<r} ^ q j m pjj eg ^ na ^ \ cannot be discrete. 
10 In particular, A cannot contain isolated points. 
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• Where does the scattering dynamics change its stability property? We 
call a bifurcation value v r y lt £ V a critical value, if the stability property 
as function of ji is not locally constant at ji = iVit, or more precisely if for 
all neighbourhoods U of v r y lt the dynamics rjE^ is (weakly or strongly) 
chaotic for some ji = /i £ U and stable for some other ji = jii £ U . 

• How can E„ be characterized at the bifurcation values v £ V and especially 
at the critical value(s) ? 

In the systems we are going to study the transitions are of the type where the 
restricted dynamics t\Yj v is stable for v < v r y lt and weakly chaotic in an open 
interval (iVit, ^crit + e)- We call such a transition to chaos abrupt or sudden if 
t\Yj v .. is stable, we call it a cascade transition if t\Yj v .. is unstable, if furthermore 

I ^crit " I ^crit " 

the bifurcation set V H ( — oo, iVit) has only isolated points with the exception of 
one accumulation point v r y lt and if Ej, crit _ e is a finite set of cycles for all e small 
enough. The notion of a cascade transition is of course motivated from similar 
transitions in the context of dissipative dynamical systems, e.g. the period- 
doubling cascade of the logistic map. 

3.3. Transition to chaos of the double horseshoe. The bifurcation set of 
the double horseshoe family can be described as follows: 

Proposition 1. The subshifts E^ = {a £ {0, if : P(V(a)) < fNi £ Z} have 
the following specification: 

(17) E^ = {a £ {0, if : (a 8 ) 8 >, < fx V/ £ Z} 

The bifurcation set V(cr, T) C [0, 2] can be characterized by a maximality property 
w.r.t. the 1-sided shift (Bernoulli-shift) c_b((/U 8 ')^o) := (/^Oi^i-' 

(18) V(a,r) = {/i£{0,lf°; /i = m°ax4( / i)} 

Proof. The first statement is trivial, (i) If v £ [1,2] has the form stated in 
the proposition, choose a = 0u. We have P(a) = v because of the maximality 
property of v. Hence any E^ with ji < v does not contain a, so that Ej, is not 
local constant at ji = v. 

(ii) If ji < 1 then symbol '1' is not allowed, S^/cr = {0}. If ji > 1 does not 
have the maximality property then there is no a £ S^/cr with P(a) = ji because 
otherwise for an a £ a we would get P {cr ks {a)) > ii- 5 /2 and a l £{ji - 8) > fi + S 
for any S > small enough, which implies P(o~ H+ks (a)) > /i. □ 
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Another characterization of the bifurcation value is given by the topological en- 
tropy of the trapped set. Define a map 7 : [1, oof— > V(cr, T), (3 1— > a with 

n-l 

(19) a n :=[(3 n (l-J2^(3- 1 )] 

where n > 1 and [x] denotes the largest integer strictly smaller than x. Then one 
has: 

Proposition 2. (1) 7 is bijective and monotone increasing. 

(2) 7 _1 is continuous (but not 7J. 

(3) log 7 - 1 (a) = ^ op ( ( 7| So ). 

(4) a 1 — y H top (a\z a ) is a Cantor function (devil's staircase). 

Proof. Part (3) is shown in [20] by using a renewal theorem. The rest can be 
concluded from the fact that the following three assertions are equivalent: 

(i) a = 7 (/3) 

(ii) For alii > 1: < a, < [/?] and EZi a ^~' = 1 and ^(a) < a for a11 
m > 0. 

(iii) for all m > 1: 1 - f3~ m < J2T=i a ^~' < 1 - 
□ 

The onset of chaos can now be described as: 

Corollary 2. The bifurcation set V(cr, T) consists of a discrete part V D [0, 1] = 
{0} and a Cantor set V D [1,2]. The bifurcation values v £ V can be classified 
according to the number of consecutive 'l's allowed in the corresponding Yj v : 
/i < 1 ji £" V; S^/cr = {0}: no symbol '1' allowed 
li = l 1 £ V; £i/<r = {0, 010}: 

first appearance of (isolated) symbol '1'; 
l</i<l.l n £ V n = 1.0 ro 10 ri . . . W* . . . 

with r = (r 8 )°^ = min{(r 8 ) 8>J : j £ N }; 
ji = l.l 1.1 £ V ; first appearance of a pair '11'; 
etc. 

/j = 1.1 = 2 1.1 £ V; Eu = {0, 1} Z : full shift. 

Proof. Vfl [1, 2] is disconnected by numbers whose binary representation contains 
a segment V with r > large enough. Its complement is obviously open, and by 
proposition 1 none of its points is isolated. □ 
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Corollary 3. The onset of chaos of the family (<r, E^) is sudden at u cr i t = 1 and 
for all fx > u cr it the topological entropy H top > and the dynamics (<r, E^) is 
strongly chaotic. 

Proof. All points in E^, ji > 1, are accumulation points, because one can exchange 
the tail of any a £ E^ by r 10 for r > large enough. The periodic sequences are 
dense because one can imitate any a £ E^ up to any indices n > 0, m < 0, then 
add a segment r 10 r , r > large enough, and repeat the whole periodically. A 
dense orbit is constructed by segments of periodic sequences which are connected 
by segments r 10 r , r > large enough. □ 

3.4. Transition to chaos of the sawshoe. The corresponding study of the 
more complicated sawshoe yields: 

Theorem 3.2 (Troll 91). The transition to chaos of the sawshoe family occurs 
via a cascade, whose bifurcation values in symbolic parameter space are 

tt^ := • , 
tt*- 1 ) := 1-^22 , 
TT {n) := 1 • -10 n -12 for n > , 

(20) ir {co) := 1 • -10 

and which accumulates in the real parameter space towards the critical line 

TT 2 

(21) d o V = {(k, 6 max ) £ V; k > 7.2 and k 6 max = — — } 

"max 

The dynamical systems on the critical line are unstable with H top = and the 
cycles are dense in the set of trapped orbits. 

To characterize the transition in terms of structural complexity we need the 
following representation of the periodic regime of the dynamics: 

Proposition 3 (Troll 91). Each cycle a £ Per(E 7r oo)/<r can be generated from 
the cycles and -11 by the transformation semi-group spanned by the following 
maps operating on E^oo ; 

(i) reflection map: 

(22) t R (a) := -a 
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(ii) basic inflations operating as replacement rules for single symbols : 



(23) 
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T+ 
J [2] 



1, 







0-11 


% ■■ 




-1 




0-11 




1 


l-> 


2 -ll 




o, 


1 


l-> 


2 -ll 


1, 





l-> 


0-11 


7 [2] : 




-1 


l-> 


0-11 




1 




0(-ll) 2 


o, 


1 




0(-ll) 2 



(iii) repetition map 

Each a £ Per(E 7r oo) ; a ^ ; can 6e formed by concatenating blocks B r [ a ] 
with fixed block class a £ {1,2} and varying block order r £ Z\{0} ; where 
B\ r \[i] := |r| -ll ; B\ r \ [2 ] := 0(-ll) |r| and B_\ r \ [a] := --B| r | [a] . 
The repetition map A operates now on the constituent blocks B r [ a ] as fol- 
lows: 



(24) A : B r[a] 

For a = we define A(0) := 0. 

The way orbits are generated from the fixed point by this semi-group is 
illustrated in figure 7. 



4. Complexity : Subshifts and Languages 

The problem addressed in this section is how to characterize the transition to 
chaos of pruned systems. The onset of chaos in such systems indicates features 
of a phase transition, seen for instance in the topological entropy as a function 
of the family parameters (cf. proposition 2 and [16]). It was conjectured (cf. 
[4]) that high level computation occurs at phase transitions. One indication 
is the appearance of infinite subshifts as a consequence of the divergence of a 
"grammatical" correlation length, which specifies how many predecessors of a 
symbol must be known to determine whether it is allowed or forbidden. It seems 
to be natural to interpret such a singularity at a critical point as indicating a 
complexity jump which requires a new and more complex production "mode" of 
the pruned SD. The theory of automata and formal languages offers concepts 
to make this more precise. Appendix C.l summarizes some basic notions and 
notations. But first we have to show how to connect the (infinite) sequences of a 
SD with the finite words of a formal language. 
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FIGURE 7. Example of the generation of a cycle by inflation. The 
codes of the cycles are given both in a 3-symbol generalized 2:1 SD 
and a 6-symbol 1:1 SD. The centers of the potentials are marked 
by cogwheel symbols. 
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4.1. Dynamically generated languages. Here we identify those languages 
which are generated by subshifts. These are the associated central language 
£(E) which consists of all non-empty segments of permitted sequences and the 
associated cyclic language £ C (E), which consists of all full periods of a cycle. 
Observe that £ C (E) contains together with any word w all its repetitions w\ 
i > and all their cyclic permutations. 

An important property of these dynamically generated languages is that £(E) 
determines E uniquely and so does £ C (E), provided E is a cyclic subshift (def. 
cf. appendix). Therefore different subshifts can still be distinguished on the level 
of their associated languages. It is possible to describe all dynamically generated 
languages on the language level. For that purpose one introduces the center and 
the cyclic center of a language: 

(i) The center C(L) of a language L consists of all biextensible words: 

(25) C(L) = {a £A + ;for all N > there are words x,y £A* 

of lengths lg(x),lg(y) > Ns.t.xay £ L} 

(ii) The cyclic center C C (L) is 

(26) C C (L) = {a £A + ; V repetitions i > V cyclic permutations 7r(a 8 ) 



We call a language L central if C(L) = L and cyclic if C C (L) = L. One sees from 
this definition easily that the associated central language of a subshift is central 
and that the associated cyclic language of a cyclic subshift is cyclic. The converse 
statement is also true, so that 

Proposition 4. There is a bijective correspondence between subshifts and central 
languages and between cyclic subshifts and cyclic languages. 

4.2. Complexity Transition. We measure the structural complexity of any 
language L by its minimal computational model, in particular by the simplest 11 
Chomsky class it belongs to. We define the Chomsky complexity x(L) of any 
language L as (-1) times 12 the type number of the simplest machine accepting 
it. A complexity transition is then roughly a change of the Chomsky complexity. 
However, we prefer to consider a particular symbolic dynamics only as a real- 
ization of a given discrete dynamical system. There are many other equivalent 

11 Observe that we chose the Chomksy type number to increase with decreasing complexity. 
12 because we prefer the reversed order; 




TRUNCATED HORSESHOES AND FORMAL LANGUAGES 



25 



symbolic dynamics. We choose a particular one among them by convenience in 
much the same way as we choose a convenient coordinate system for a geometric 
problem. Therefore, we have to make sure that a local complexity change occurs 
in each equivalent family of symbolic dynamics. Indeed, this independence of the 
symbolic dynamics is a general theorem (cf. [17, 19]). 

Definition 4.1. The pruned system (r^, E |U ) |U£ j shows a complexity transition at 
the critical value /i cr ;t £ J if for any equivalent symbolic dynamics (<r, E |U ) |U£ j, the 
Chomsky type function n \— > x(jC(E^)) is not locally constant at ji = /i CT it, i-e. 
assumes more than one value in any neighborhood of /i cr it- 

Theorem 4.1 (Troll 92). Any semi-conjugacy $ between two subshifts does not 
increase the Chomsky complexity of the associated (cyclic or central) language. 
In particular, if $ is even a conjugacy then: 



4.3. Complexity transitions in truncated horseshoes. We state the main 
results for our horseshoe models. First for the double horseshoe: 

Theorem 4.2 (Troll 92). Suppose v £ Vf~l]l, 2] . Then both the associated cen- 
tral language jC(E^) and the cyclic language £ C (E |U ) are regular languages if and 
only iff G]l,2] is rational. The associated languages are finite complement lan- 
guages (i.e. E^ is a finite subshift) iff v has a periodic binary representation 
without leading transient part. 

Corollary 4. A complexity transition occurs for the truncated double horseshoe 
simultaneously with the onset of chaos at u cr i t = 1. 

For the sawshoe one can state: 

Theorem 4.3 (Troll 91). The parameter family of the associated formal (ei- 
ther central or cyclic) languages has the following properties: 

2 

(a) in Vo = {(£;, b max ) £ V; k > 7.2 and kb max < — | — } the language is 
finite (hence regular); 

(b) on d Vo the language is not context free but context sensitive. 



(27) 



x(£(S A )) = x(£(S fi )) 






assumed for the second equality that E^ is cyclic. 



ae 
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Appendix A. Basic definitions in symbolic dynamics 

We discuss here dynamical systems whose dynamical map operates on se- 
quences. We permit sequences over some finite set A of abstract symbols with 
cardinality $A > 1 and define the full shift over A as the complete set of sequences 

We can make a metric space out ofA z by defining distances between sequences: 

(29) ^M) = E%r 

where a, b gA z , S xy = if x = y, 1 otherwise. 

Among the continuous maps operating onA z shift maps and cellular automata 
(CA) are of special interest because they are compatible with the group structure 
of the index set Z: 

Denote by a% the (right) shift map onA z : (<r(a)) 4 - = a 8 ' +i . Usually we write just 
a for (T%. A subshift is a closed cr-invariant set inA z . E is a finite subshift if it is 
determined by a finite set F of forbidden symbol strings, which are not permitted 
to appear in sequences. A subshift is Markov if the immediate predecessor deter- 
mines whether a symbol is allowed or forbidden. Such a Markov rule is usually 
given in the form of a transition matrix indexed by the alphabet and showing an 
entry 1 wherever a pair of symbols is allowed and otherwise. We call E cyclic 
if the periodic (w.r.t. a) sequences are dense. A dynamical system given by s 
subshift E cA z and the shift map <r|E — > E is called a symbolic dynamics. 

A natural generalization leads to cellular automata (CA). A (1-dimensional) 
CA is any continuous map r onA z which commutes with the shift map. 13 

A way to compare dynamical systems topologically is offered by semi-conju- 
yacies. We define them here just for shift maps as any continuous map $ between 
subshifts E^, E^ which intertwines between the respective shift maps <r^, cry. 

(30) $ : E A E fe , ay o $ = $ o a% 

If moreover, $ is a homeomorphism i.e. bijective and also $ _1 continuous, then it 
is called a conjuyacy (between shift maps). It may be interpreted as a continuous 
change of variables which identifies the shift maps on E^ and E^. For instance, 

13 Since any CA is automatically uniformly continuous an equivalent description of a CA 
(cf. [8, 9]) can be given by a local function f T : A 2R+1 A with the property (r(s)) 8 - = 
/ t (s 8 __r, . . .,s i+R ) for any s 
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finiteness and the Markov property of subshifts are equivalent in this sense, i.e a 
subshift E is finite iff it is topologically conjugate to a Markov subshift. 

Appendix B. The inverse scattering problem for the sawtooth 

DEFLECTION 

B.l. Statement of the inverse problem. We are looking for an attractive 
central potential V : ]0, oof— > IR whose deflection function By on ]0, & ma x[U]& max , oo[ 
(allowing for singularities at and 6 max ) for a given asymptotic particle energy 
Eq > is the sawtooth k > 0. We require the following properties: 

(1) V <0,V>0 

(2) V(r) = if r > 6 max , 

(3) V continuous, 

(4) V smooth (C 2 suffices) in ]0,6 max [. 

The forward problem of finding the deflection function By of a given (not necessar- 
ily attractive) central potential V of range 6 max and a given asymptotic particle 
energy E > can be easily solved from the Hamilton equations to yield for 
impact parameter < b < 6 max the possibly infinite value: 

POO A r 

( 31 ) 7 \ , 

v ' b f b ^ dr 

= —2 arccosf- ) + 2 / , = 

o max J rrain (b) r 2 ^J(i _ y^l)b- 2 - r- 2 



where r m ; n (6) is the classical turning point (minimal r value) of the trajectory. If 
V'(r) = (zero force) can be excluded, r m ; n (6) coincides with the most remote 
point satisfying r = 0, i.e. it is determined as the largest zero of the radicand 
(l-^b-i-r- 2 : 

(32) r min {b) = max{r > 0; 1 - V(r)/E = b 2 /r 2 } 
For b < we set B v (-b) := -B v (b). 

B.2. Solution of the inverse problem. In order to present the solution for- 
mula conveniently we define a function p: 

(33) p : [0,6 max ] -> [0,6 max ], s _^ s e"^ v 7 ^- 2 
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As one easily verifies (see proof below), it is strictly monotone increasing, surjec- 
tive, continuously differentiable and it has no critical points in [0, 6 max [- Hence, 
p has a continuously differentiable inverse function £. 



Theorem B.l. A solution of the inverse scattering problem B.f is given by: 
(34) 

V Zk : [0,oo[-> E ", r ^ 



E (l - e 2fc /^V 6 max-C 2 (r)^ Z J o < r < 6 n 
ifr> 6 max 



The solution is unique in the class of potentials satisfying V(r) > 1 — 6^ ax /r 2 Vr £ 



]0,6 max [- 



A parameter representation of this potential is given by 

(35) {(s e -k/WZ™->\ i?o(l-e 2fc/ ^ 6 — " s2 )) ; s £ [0,6 max ]} 

In the proof we shall proceed as follows. First we shall show that equation (34) 
solves the inverse problem B.l. Afterwards, we shall prove uniqueness. 

B.2.1. Existence. 

Proof. That p is strictly monotone increasing and has no critical points in [0, 6 max [ 
can be seen by differentiating equation (33): 



/q«\ // \ ks 2 + TT\/b max 2 - s 2 

(36) p(s) 



,k/7r y/b^S-s 2 tt Jb 2 - S 2 



Hence p and its inverse £ are admissible integral substitutions in equation (31). 
Indeed they are even analytic functions, so that Vz K is also analytic in ]0,6 max [ 
and, obviously, continuous in [0, oof. Substituting s = £(r) and (34) in (31) yields 

(37) 9(b) = kb 

for < b < 6 max . □ 
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B.2.2. Uniqueness. 

Proof. Let V be a solution of the inverse problem B.l. We first observe that By 
is most naturally a function of b~ 2 . Therefore we introduce 

(38) y:=b-\ y :=b 2 max , 0(y) := v (b), and v(r) : = 1 V(r) 



Ei 



o 

Next we make a substitution in the integral of equation (31), whose purpose is to 
remove the implicit dependence of the integral limits (i.e. r m ; n ) on the potential 
V: 

(39) x(r) := 

r z v{r) 

Lemma 2. If v(r) < b 2 nax /r 2 y r G]0,6 max [ ; then the variable transformation x is 
a positive and monotone function bounded by r~ 2 . 

Proof. The critical points of this transformation satisfy: 

... 2r v(r) + r 2 v'(r ) 

40 x'{r) = V [ , , . y J = 

v ' v ' r A v 2 (r) 

i.e. 2r v(r) + r 2 v\r) = 0. The existence of a critical point r is obviously equiva- 
lent to the existence of a constant C > 0s.t.u(r o ) = C/tq and v'(r ) = — 2C/?"q, 
i.e. v is locally inverse quadratic at any critical point. In order to show that 
equation (39) defines an a monotone function it suffices to show that x does not 
have a critical point in the open interval ]0, & max [- 

Lemma 3. Let < b < 6 max be the impact parameter of a trajectory with clas- 
sical turniny point < r < b mSLX , where the potential v is assumed to be locally 
inverse quadratic. Then the deflection function of v has a singularity at b and 
is consequently unbounded. 

Proof. One way to see this is to set the Taylor expansion of v around r into 
(31). One finds a nonintegrable singularity (r — r ) _1 . A more physical argument 
is given by the observation that criticality at r is equivalent to the equality of 
the potential's gradient at r and the centrifugal force for curvature radius r . 
Together with the turning point condition r(r ) = and the uniqueness of the 
initial value problem of motion in a potential field, this yields a trajectory which 
ends as circular motion, winding infinitely often around the potential's center. □ 
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Now, consider the turning point function 

(41) r min : [0, 6 max ] [0, 6 max ], b ^ max{r; r 2 u(r) - 6 2 = 0} 

Although r m ; n (0) = and r mm (6 max ) = 6 max , not every critical point r need 
to appear as a turning point for some impact parameter 6, because the turning 
point function need not be continuous. However, since v was supposed to be C 2 
in ]0, fe max [, the only way a discontinuity may arise is via a tangent bifurcation at 
a local extremum of 

(42) v :]0,oo[-» [0,oo[, r ^ r 2 v(r) 

Two cases of discontinuities of the turning point function r m ; n can be distin- 
guished: 

case 1: There is a critical point at which v assumes a value < & max . But 
then there must be relative minima among the critical points. Since v is 
C 2 there is a lowest minimum. By construction it is both a critical point 
and a turning point. Thus the conditions of lemma 3 are satisfied, so 
that trajectories with infinite total curvature occur, in other words the 
deflection function must be unbounded. However, the sawtooth deflection 
function is bounded. Hence, in this case the integral substitution x of 
(39) is a monotone function and consequently admissible. 

borderline case: Suppose there is a critical point of v on the line & max 
that is the accumulation point of a sequence of turning points. This would 
again lead to a limit circle, so that this case is also excluded because By 
would be unbounded. 

case 2: v assumes values > 6 max and all critical points of v in ]0,6 max [ 
assume values > & max . In this case there is a < r* < b mSLX s.t.r 2 v(r) > 
b 2 Vfc > r > r* 

"max v "max ^ — 

Positivity of x is immediate, because we assumed V to be nonpositive, so that 
v(r) > 1. □ 

Introducing the function 

(43) gix) = -y^±^l 

r ax 

we get 



(44) §(y) = -2 arccos(^) + 2 f dx 9[X) 
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The important point is that y = b~ 2 is both a limit of integration in the expression 
for the deflection function and a natural choice of argument for the deflection 
function itself . This makes it possible to write our inverse problem as an Abel 
integral equation. We first write the sawtooth as a function of y : 

(45) z k (y) = 4 

so that the integral with Abel kernel (y — x)~ l l 2 and upper limit y equals a 
function 6 of y: 



(4o) 0{y) := h arccos(- — J = I ax 



2^y y/y Jyo Vv ~ x 

This is an Abel integral equation for g which can be solved by the following 
inversion formula, because 6 is smooth in ]0,6 max [: 

(47) 9 (x) = ^ I -feL + r d§{y) dy 




2 Jy dy ( x _ y )2 



{x-y )' »>y [x-y 

This yields 



, , 1 k k Jx — yn 
(48) g( x ) = 1 -—^ — 

2y/x 2 7T y/x - yo yftjo 2 TT X y/t^ 

Now, in order to obtain v one uses equations (39) and (43) to find the following 
expression for w(x) := v(r(x))\ 

dln(w(x)) 1 2g(x) 



dx 

( 49 ) i, r- 

kJyo 



ttx2 y/x - y 



This equation has to be integrated. We use w(y ) = v(b msix ) = 1 (because 
x{b max ) = y ) to get: 

(50) lnWx)))- 2fcv ^" 



^ v x vyo 
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Hence 

(51) w(x) = e^/Wvo 1 -- 1 

Now, we are ready to take the last step in calculating the wanted potential. 
Observe that 

(52) r(x) = 1 = ,-1/2,-^V!/,- 1 -- 1 



is the inverse of x(r). Although it is not possible to solve equation (52) explicitly 
for x, we can easily give a parameterization of the potential curve (r, Vz K (r)) } 
r G supp(V), in terms of x G [j/o, oo[: 

(r(x), Vz K (r(x))) = (r(x), E {1 - w(x))) 

(53) = (x-^e-V*^ 1 -*- 1 , E (l - e^V^ 1 -*" 1 )) 

Monotone reparameterization by s = x~ l l 2 yields the assertion. □ 
Figure 3 shows this potential, which has the following properties: 

B.3. Properties of the potentials Vz K - Each Vz K is continuous at 6 max , be- 
cause w(y ) = 1 (see above). At r = 0, which corresponds to x — > oo, one 
finds: 

(54) y min = V Zk (0) = E (l- e 2kb -^) 
The derivative of Vz K is given by 



(55) {r(s),V ZK (r(s))) = lr(s),E 



2 g3 k/lT -^/femax 2 -S 2 J~ 



k s 2 + 7T yb, 



One sees from this expression that the potential Vz K has independently of k and 
Eq a unique critical point {V z , = 0) at the center r = s = 0. At the border of 
the potential's support (r = s = 6 max ) the left derivative is also independent of 
k, namely V z ,(6 max ) = 2£ , /& m ax ^ 0. This is necessary to get a non vanishing 
(indeed the maximal) deflection just at the border of the potential. 
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The influence of a variation of the parameters k and 6 max on the shape of the 
potential Vz K can be studied through the second derivative of Vz K - 

(56) (r(s)^ K (r(s))) 

1 2e 4fc /^ m ax 2 -s 2 k fb max 2 7 r 2 -3k 2 s 4 -4kTrs 2 Jb mi J-s 2 ) 
r(s),E V J 



y yks 2 + tt \Jb max 2 - S 2 j j 

Ar r = one finds V£ K (0) = 2kE /{Trb max ) exp(4A;6 max /^) ^ 0. The zeros of V% K 
are determined by the equation 

(57) 7r 2 b 2 max -3k 2 s*-4k 7 rs 2 (b 2 max -s 2 y/ 2 = 

The left hand side of this equation defines a function f(s) with /(0) = 7r 2 &max > 0- 
If k > is too small then f(s) > for all < s < 6 max < 1/2 and does not 
have any zeros. If however f(b max ) = (?r 2 — 3A; 2 6^ ax )6^ ax < 0, i.e. if 

TT 

(58) A; 6 max > -j= 

then V^-' , has zeros in the interval [0,6 max ]. For the third derivative one finds 
again a zero at r = 0: V^' r (0) = 0. The graphs of V z , for some typical A;- values 
are shown in figure 8. 

B.4. The unrestricted solution space of the inverse problem. If one drops 
the condition that v must not have only critical points above 6^ ax then the so- 
lution of the inverse problem is no longer unique. From the former sections we 
know that any solution cannot have a critical point below 6^ ax , because otherwise 
the deflection function would be unbounded. Let r* = inf{0 < r < 6 max ; v(r) = 
^m ax }- Then we can perform the integral transformation (39) in the interval 
[0, r*\. The value r* and the potential in [r*, 6 max ] can be chosen freely > b 2 m ^Jr 2 
but in the neighbourhood of r* where we require smooth matching and at 6 max 
where we require continuity. Going through the calculations above one finds that 
all solutions have the general form : 

b 



(59) w^ n (x) = ^^) e 2 /<o^* ( < ( ^ OBh V^ 
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V /E 



k = 4, 8, 12, 20 




Figure 8 



First derivative of the potential Vz K for some values of k. 
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where 
(60) 



R r *{y) 



dr 



r\yv{r) 



B.5. Motion in the potential Vz K - Integrating the equations of motion for 
impact parameter < b < 6 max yields the following trajectories in polar coor- 
dinates (r, (f>) as functions of the parameter s £ [6, oof (there is a singularity at 
b = b m ^)\ 



<f>(r(s)) 



(61) 



r(s) 



dz 



mm O 

z l 



arccos 



arccos 



arccos 



arccos 



U^M^_ 2 



En 



it arctan(2 

Ztt v 



-si), 



6 2 +6^ ax -2 S 2 

if b < s < s 



sing 



\p2b 



b_k 

4 



if, 



^sing 



it arctan(2 

Ztt v 



-si), 



2 



if ^sing < S < b n 



if s > b n 



b_k 

2 



where s s ; ng = (b 2 + 6 niax ) 1/ ' 2 /v / 2 is the singularity in the argument of the arctan. 
The trajectories are continuously differentiable everywhere, e.g. for r = s = 

brnav- 



(62) 



lim 4>'(r) = lim </>'(r) 

r->6 m ax+0 r->6 max -0 



b m ^Jb 2 max - b 2 



Some sample trajectories are shown in figure 4. ae 



Appendix C. Formal Languages and Turing machines 

In this section we introduce notions which will serve to characterize the struc- 
tural complexity of subshifts. 
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C.l. Basic definitions. We start from a finite setA of abstract symbols, which 
will be called alphabet in the following. Define the set of all words or (finite) 
stringsA* = \J™ =0 A n . ^ 1S called the total language overA and together with the 
concatenation operation becomes a semi-group. Its unit element is the empty 
string $. Powers correspond to repetitions of symbols, the length lg(s) of a string 
s gA* yields a formal logarithm. Any subset L of the total languageA* is called 
a language, e.g. A + := A* \ {$}. A word x £ L is a subword or segment of a word 
y £ L, written x -< y if x appears within y, i.e. if there are a, b £A*s.t . y = axo. 

C.2. Automata. Automata are mathematical models of devices that process 
information by giving responses to inputs. Formal language theory views them 
as scanning devices or acceptors able to recognize, whether a word belongs to a 
given formal language. There is a hierarchy of automata corresponding to the 
complexity of their structure, in particular of their memory. Let Z be a finite 
nonempty set of internal states Z{ } among them an initial state z 0} and a nonempty 
set of final states Zf C Z. Let A and VV dA be two (as usual nonempty, finite) 
alphabets, called the input alphabet and the working alphabet respectively. The 
working alphabet contains at least one extra symbol, the blank. We list the 
following basic types z = 0,1,2,3: 

(0) Turing machine (TM): A (deterministic) TM given by the tuple 
M = (Z,A, VV, u, r, /i, zo, Zf) is defined by three transition maps, an over- 
writing map u : Z X °W — > \V, which overwrites the symbol being read, 
a next state map r : Z x\V — > Z , which defines the machine's transition 
to a new internal state, and a head moving map fi : Z X °W — > {L, i?}, 
which moves the reading head of the machine by one position to the right 
or left. The language accepted by M is the set of words w £ A* which 
make M, starting from the first letter of w, reach a final state. 

(1) Linear bounded automaton (LB A): A LBA is a Turing machine 
whose work space, where its head is allowed to move, is restricted to (or 
equivalently bounded by a linear function of ) the space the input word 
uses on the tape. 

(2) Pushdown automaton (PA): A PA uses a stack (first in - last out 
memory) instead of a tape as memory. 

(3) Finite deterministic automaton (FDA): A FDA does not have 
an extra memory. It is given by a single function, namely the next state 
or transition function r, i.e. the FDA performs the following actions: It 
scans the letter s 8 - of the word w, while in an internal state z £ Z, it 
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performs a transition to the new state z' = t(z } and then it moves to 
the next letter This is repeated until the last letter of w is reached. 

If the last state is a final state, the word is accepted. 

A hierarchy of automata corresponds to the Chomsky hierarchy of languages, 
namely the recursively enumerable (0), context sensitive (1), context free (2) 
and regular (3) languages, in the sense that any class i language is recognized 
by a class i automaton and conversely any class i automaton produces a class i 
language as output. 
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